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In this talk I review the current status of jet formation in direct numerical simulations of black-hole 
accretion disks and magnetospheres. I address the following critical questions: What constitutes 
the jet? What is the launching mechanism? Where is the launching point of the jet? What is 
the Lorentz factor? What is the opening angle? How is the jet collimated? Just as importantly, I 
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Let me begin by saying a few words about what I will not cover in this review. I plan to stick 
to a narrow interpretation of the title given to me by the Scientific Organizing Committee. Firstly, 
I will be dealing mostly with direct numerical simulations; I will not present any observational 
results and will only briefly mention one or two well-known theories. Secondly, I will specifically 
focus on numerical simulations of the formation of jets, ignoring the large body of literature on 
the evolution and stability of jets and their interactions with their environments, where the jets are 
taken as an input condition. I will instead carefully focus on simulations that include complete 
treatments of the accretion disks and the region right around the event horizon of the black hole, 
where the jets are actually launched. 

A convenient place to start this review is to look at where this field was at the time of the 6th 
Microquasar Workshop in 2006. That also happened to be the year that two of the main papers on 
numerical simulations of black-hole jet formation were published, one by McKinney [13] and the 
other by Hawley [8]. After summarizing their work, I will discuss some of the important things 
that have been learned about jet formation in MHD simulations in the two years since then. I will 
finish by discussing some possible future directions for this research, some of which is already 
under way. 

1. What Constitutes the Jet? 

Let me first be clear about what is meant by the term "jets" when referring to numerical 
simulations, as it is not obvious at this point how they are related to the jets observed in Nature. 
In the simulations of McKinney [12] and Hawley [8], the jets are identified as unbound (ut < — 1), 
outflowing (u^ > 0) material, where is the 4-velocity of the fluid. It is actually fairly remarkable 
how similar the results of such numerical simulations are, despite a number of differences in the 
particulars of the codes that are used. A quick comparison shows that the numerical simulations of 
both groups produce qualitatively and quantitatively similar final states, which include prominent 
jet-like features. The various features identified in these simulations are illustrated in Figure 1. 

For most of the rest of this review I will distinguish between the two different jet components 
shown in Figure 1: the Poynting-flux jet and the matter-dominated or funnel-wall jet. Despite 
the strong similarities in their results, McKinney and Hawley chose to focus on different aspects 
of their results. Most of McKinney's analysis is concentrated on the Poynting-flux jet, whereas 
Hawley generally discusses the matter-dominated jet in more detail. This will make it easy for me 
to keep the two separated in this review. 

2. The Launching Mechanism 
2.1 Poynting-Flux Jet 

Figure 2 shows magnetic field lines overlaid on a pseudocolor plot of density from the sim- 
ulations of McKinney [13]. The important things to note about that figure are that, although the 
magnetic field lines appear tangled and chaotic inside the disk, as expected from the action of the 
magneto-rotational instability (MRI), inside the nearly evacuated funnel region the field lines are 
well ordered all the way down to the event horizon of the black hole. In fact, the magnetic field lines 
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Figure 1: Left panel: Schematic illustration of the main dynamical features seen in the accretion disk 
simulations of Hawley [8]. The jet is defined as those regions that are unbound and have positive radial 
momentum. Right panel: A more quantitative delineation of the different regions from the simulations of 
McKinney [12]. The material inside the green contour, corresponding to the main disk body, is weakly 
magnetized (j3 = Pgas/Pmag ^3). The material inside the violet contour, corresponding to a magnetized 
corona, is near equipartition (1 < j3 < 3). All of the material inside the black contour is unbound {ut < —1) 
and would be called "jet" material. Inside the red contour, the magnetic field energy density dominates even 
over the rest mass density of the fluid {B^ / p > 1). 



have the characteristic appearance of the split-monopole configuration described by Blandford & 
Znajek [3], suggesting a possible launching mechanism. 

It is worth emphasizing at this point that the initial configuration for McKinney's simulation 
did not include any large-scale magnetic fields. Instead the simulation started with only a weak 
dipole magnetic field in the form of poloidal loops contained entirely within the disk (see Figure 
8 below for an illustration). The fields that subsequently fill the evacuated funnel and power the 
Poynting-flux jet arise solely from these initial seed fields and the evolution of the disk. This is 
important because it means that powerful jets do not require large-scale, ordered magnetic fields 
be fed in from large radii through the disk. 

Although the field configuration in Figure 2 is strongly suggestive that the Poynting-flux jets 
are powered by the Blandford-Znajek process, it is not conclusive. Fortunately the case is made 
much stronger by the careful quantitative comparisons McKinney & Gammie made between their 
simulations and the predictions of the Blandford-Znajek model [12]. Examples of the results of 
such comparisons are shown in Figure 3. Because the Blandford-Znajek model does not include a 
realistic accretion disk, one must be careful when comparing it to numerical simulations of disks, 
such as the ones being considered here. It is only realistic to expect the Blandford-Znajek solution 
to apply in regions where the magnetic field is dominant, as in the nearly evacuated funnel region. 
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Figure 2: Magnetic field {black lines) overlaid on a pseudocolor plot of the log of density from the simula- 
tions of McKinney [13] 
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Figure 3: Left panel: Square of radial field on the horizon from the time- and hemisphere-averaged simula- 
tion data of [12] {solid line). The dotted line shows the Blandford-Znajek [3] perturbed monopole solution 
with the field strength normalized to the numerical solution at the pole. Right panel: Electromagnetic energy 
flux on the horizon from the time- and hemisphere-averaged data of [12] {solid line). The dotted line shows 
Blandford-Znajek' s spun-up monopole solution. 



From Figure 3 we see that in this region (0 < 7r/4) the numerical simulation agrees quite well 
with the model. This would seem to confirm that the Poynting-flux jets are indeed powered by the 
Blandford-Znajek mechanism (magnetic field lines threading the ergosphere of the black hole tap 
the rotational energy of the black hole and transport that energy to large radii). 

2.2 IMatter-Dominated Jet 

The matter-dominated jet has a much different launching mechanism. In some ways it is even 
easier to understand, particularly if one can understand Figure 4, taken from the simulations of 
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Figure 4: Color contours of the azimuthal average of the total pressure (gas plus magnetic). Radial momen- 
tum flux in the outflow is shown in white contours. Both sets of contours are linear scales and come from 
the simulations of Hawley & Krolik [8]. 



Hawley & Krolik [8]. Notice in the figure that there is a strong pressure gradient along the coronal 
wall that points oblique to the momentum-flux contours. Breaking this pressure gradient into a 
component perpendicular to the momentum-flux contours and a component parallel, we can easily 
see that the pressure gradient serves two roles: 1) the perpendicular component pushes matter 
into the jet (mass loading); and 2) the parallel component accelerates matter along the jet. These 
strong pressure gradients line the entire funnel wall, so there is a very extended mass loading and 
acceleration region. 



3. The Launching Point 

We can confirm that the acceleration of the matter-dominated jet continues all along the funnel 
wall by looking at the plot of mass outflow rate shown in color in Figure 5. Note that the mass 
outflow rate continues to increase out to the outer radial extent of this figure (r ^ 15rG). Figure 5 
also shows the equipotential surfaces {black lines) of the flow. Motion of gas away from the black 
hole across the equipotential lines inside r < 6rG, means that the jet material is climbing out of the 
potential well of the black hole. As shown in Figure 4, the driving force behind this acceleration 
is the strong gradient in total pressure near the hole. Another important point illustrated in Figure 
5 is that outside of r > 6rG, the matter-dominated jet flows parallel to the equipotential surfaces, 
meaning that the acceleration must not be attributable to centrifugal forces. This strongly rules out 
a launching mechanism such as the Blandford-Payne wind model [2]. 



4. Lorentz Factor 



Now we turn to consider characteristics of the jets in these simulations that might be compared 
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Figure 5: Mass outflow and effective potential contours in simulation KDR from [8]. Here, the color 
contours are the logarithm of the mass outflow rate, while the line contours are linear in the potential. 

directly with observations, such as the Lorentz factor and opening angle. For the matter-dominated 
jet, measuring the Lorentz factor from numerical simulations is a rather straightforward operation, 
and a range of values are quoted in the literature. Typical values are F ^ 1.4 — 1.5 [12, 9]. 

However, for the Poynting-flux dominated jet, measuring a Lorentz factor from numerical sim- 
ulations like those of McKinney or Hawley is a difficult proposition. This is because the evacuated 
funnel region occupied by the Poynting-flux jet can not be handled properly with the numerical 
techniques that are employed. This is because the evacuated funnel should really be evacuated 
(void of any matter) near what is called the stagnation point. The stagnation point, quite simply, 
is a point along each field line, inside of which matter is falling toward the black hole and outside 
of which the matter is moving off to infinity. In ideal MHD (as treated in the numerical codes of 
Hawley and McKinney), matter is not allowed to move across field lines, so it is impossible to 
replenish matter as it is drained away from the stagnation point. Eventually, the density there must 
go to zero. However, in the numerical codes used in the simulations Fve been describing, the gas 
density can never be allowed to go to zero without causing the code to crash. Therefore, matter 
must artificially be added to the simulation in the vicinity of the stagnation point and anywhere 
else the density falls below some threshold. Because the mass content of the jet is a main factor 
in determining the Lorentz factor, this artificial mass-loading can have a profound impact on the 
Lorentz factor measured in the numerical simulations. A realistic Lorentz factor could only be 
obtained if there were a valid model for the mass loading of the jet. 

5. Opening Angle 

Along with being relatively slow, the matter-dominated jets also have fairly large opening 
angles, as shown in Figure 6 from Hawley & Krolik [8]. The initial opening angle at the base of the 
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Figure 6: Opening angle of the jet in radians, as defined by the angular location of the boundary between 
bound and unbound outflow, from the simulations of Hawley & Krolik [8]. Four models are shown: the solid 
line is a/M = — 0.9, the dotted line is f2/M = 0.5, the dashed line is a/M = 0.9, and the dot-dashed line is 
a/M = 0.99. 

jet is in the range 50° ^ ^ 60°. As the matter-dominated jet moves outward, some coUimation is 
provided by the coronal walls, so that at r ^ 40rG the opening angle is down to 40°. Still, this seems 
rather too large to agree with most observed jets. However, it may be that there are mechanisms, 
such as the internal hoop stress within the jet or pressure gradients from the external environment, 
that could further collimate them at even larger radii. The final answer on the collimation of the 
matter-dominated jets will have to await simulations that extend to larger radii and run for longer 
times. 

One beneficial attribute of the simulations of Mc Kinney [13] is that they covered a much larger 
spatial extent than previous simulations of this type. The outer boundary of his simulations was at 
^max ^ 10"^ rG (as opposed to r^ax ^ lOOrG in earlier simulations), so he was able to follow the evolu- 
tion of the jets over nearly 4 orders of magnitude in spatial extent. Again his analysis concentrated 
on the Poynting-flux dominated jet, which he found starts with only a slightly smaller opening 
angle than the matter-dominated jet (0 ^ 40°). Like the matter-dominated jet, the Poynting-flux 
jet is initially confined by the coronal wall. After extending beyond the range of the corona, the 
Poynting-flux jet is subsequently confined by the surrounding matter-dominated jet. Eventually 
the Poynting-flux jet is able to collimate itself through the action of internal hoop stresses from 
the toroidal component of the magnetic field that is carried along in the jet. The self-collimation 
continues until about r ~ 500rG, when a final opening angle of ~ 7° is achieved. The gradual 
collimation of the jet is illustrated in Figure 7. 

6. Caveats 

What I've reviewed up to this point was the state-of-the-art in MHD simulations of jet for- 
mation at the time of the 6th Microquasar meeting. I now shift to consider some things we have 
learned about such simulations since that time. 
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Figure 7: Jet opening angle as a function of radius in the simulations of McKinney [13]. The horizontal axis 
is the logarithmically- spaced radius from Ivq to lO^rG- 



6.1 Initial Conditions 

Ideally one would hope that the results of the simulations I have been discussing would be in- 
sensitive to their initial conditions, particularly as the simulations are run out to longer times. Since 
these simulations start from highly artificial initial conditions (a simple gas torus seeded with a 
simple magnetic field topology) as illustrated in Figure 8, it is particularly important to look for the 
development of a long-term, steady state disk-jet solution. It is this disk-jet solution that one hopes 
is independent of the initial conditions. However, because magnetic field in Nature are divergence- 
free, they can never truly "forget" their initial configurations. In analyzing simulations of the type 
I have been reviewing, Beckwith et al. [1] discovered that, although the evolution of the disk is 
essentially independent of the initial conditions, this can not be said for the jets. By comparing 
the evolution of initially identical tori seeded with different magnetic field configurations (such as 
dipole, quadrupole, multiple loops, and toroidal as shown in Figure 8), Beckwith et al. discovered 
that a critical component in the formation of strong jets is that the magnetic field have a net poloidal 
flux. Without this, any jets that form appear to be orders of magnitude weaker, as illustrated using 
a variety of measures in Figure 9. 



7. Future Directions 

I'd like to finish this review by highlighting a couple areas where new contributions are being 
made in the area of jet formation in MHD simulations. I'm particularly interested in these because 
they are the current focus of my own research. 

7.1 Tilted Disks and Jets 

A fundamental question regarding jet formation that remains unanswered is: To what charac- 
teristics of an accretion disk system can the orientation of a jet be ascribed? Fundamentally there 
are three possible answers to this question. The first possibility is that the jet orientation is depen- 
dent on the angular momentum or spin of the black hole [3]. This would suggest that the jet should 
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Figure 8: Initial configurations of dipole {left panel), quadrupole {center panel) and multiple loop {right 
panel) field topologies. The torus for the multiple loop topology is shown slightly zoomed in to illustrate 
better the initial field structure. White contours denote magnetic field lines, color contours the gas j3 param- 
eter. Solid and dashed lines indicate opposite field polarities. 
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Figure 9: Time average of data from unbound outflows as a function of radius for models with a dipole 
{black lines), quadrupole {blue lines), and toroidal {magenta lines) initial magnetic field. These images are 
taken from [1]. Shown are mass outflow rate M {top-left panel), \ \b\\^ {top-right panel); the electromagnetic 
energy flux {bottom-left panel); and the angular momentum flux {bottom-right panel). Dashed lines show 
±1 standard deviation from the average. 
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always be perpendicular to the black-hole symmetry plane. The second possibility is that this ori- 
entation is dependent on the angular momentum of the accretion disk [10], thus orienting the jet 
perpendicular to its plane of symmetry. A final possibility is that a magnetic field, generated by or 
near the central object or carried in through the disk, dictates the jet orientation [2]. In theoretical 
work it is often assumed that all of these components are aligned, even though there is certainly 
no requirement that this be the case in real systems. Such an assumption makes it difficult to de- 
termine what is setting the orientation of the jet. From an observational perspective, the angular 
momentum vector of the central black hole is not directly observable. In fact, its orientation has 
traditionally been inferred from the orientation of the jet and its magnitude from the properties of 
the disk. This obviously creates a dilemma if one wants to observationally constrain what is setting 
the orientation of a jet. 

However, direct MHD numerical simulations of tilted black-hole accretion disks, such as those 
presented in [5], allow one to explore situations in which the black-hole spin, disk angular mo- 
mentum, and magnetic field are not all aligned. One then need only observe the orientation of any 
resulting jets in these simulations to address the question of what sets their orientation. The disks in 
tilted simulations have also been shown to precess [5] ^ so another interesting question is whether 
or not jets can be forced to precess by such a disk. 

7.2 Beyond Ideal MHD 

Perhaps the greatest shortcoming of the numerical simulations of black-hole accretion disks 
and jets carried out up until this past year, has been the very incomplete treatment of thermodynam- 
ics. This is also the area in which the simulations of McKinney and Hawley differ the most. Hawley 
and his collaborators use a formulation of the MHD equations that evolves the internal (not total) 
energy of the gas [4]. Because their code does not conserve total energy and does not include any 
terms to capture dissipative heating processes, any heat created by disk turbulence is simply lost 
from the simulation. In a sense, the entire flow is artificially cooled, although heating from com- 
pression and shocks is still treated. McKinney, on the other hand, uses a total-energy conserving 
formulation [7], so that any kinetic and magnetic energy dissipated within the disk is automatically 
recaptured as heat. However, McKinney's code does not include any radiative cooling processes. 
As such, it strictly only applies to radiatively inefficient flows such as the disk around Sgr A*. 
Recently, Fragile & Meier [6] performed the first global numerical simulations of black-hole ac- 
cretion disks that included both a total-energy formulation to capture dissipative heating processes 
and a radiative cooling function that approximated cooling due to bremsstrahlung, synchrotron, and 
inverse-Compton processes. Figures 10 and 11 show how these different numerical approaches af- 
fect the final state of the disk. The key result is that the radiatively-cooled disk is much cooler and 
much closer to equipartition than the other models. The approach toward equipartition may signal 
a transition to a magnetically-dominated accretion flow (MDAF) as proposed by Meier [11]. 

The other important point about including real radiative cooling processes in numerical simu- 
lations is that it opens the door to the possibility of making real connections between simulations 
and observations. This is an exciting new possibility! 

^Movies of the tilted disk simulations of Fragile can be seen at http://fragilep.people.cofc.edu/research/tilted.html . 
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Figure 10: Pseudo-color plots of the logarithm of temperature with contours of the logarithm of density. The 
upper-left panel is from an internal energy model (like De Villiers & Hawley); the upper-right panel is from 
an internal energy plus cooling model; the lower-left panel is from a total energy model (like McKinney & 
Gammie); and the lower-right panel is from a total energy plus cooling model from Fragile & Meier [6]. 



11 



Jet Formation in MHD Simulations 



P. Chris Fragile 




40 80 

r/r^ 



Figure 11: Main disk properties plotted as a function of radius for the internal-energy model I, the internal- 
energy plus cooling model I+C, the total-energy model T, and the total-energy plus cooling model T+C from 
Fragile & Meier [6]. Plotted variables are the density p, gas pressure P, gas temperature T, temperature 
scale-height H/r, stress parameter a, and the equipartition parameter j3. 
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8. Conclusions 

The main take-away points from this review are 

• Jets are consistently formed in MHD simulations. They do not require pre-existing large- 
scale magnetic fields. They do, however, require a net poloidal flux. 

• The jets in these numerical simulations have a 2-component, sheath configuration, with a 
slower, matter-dominated jet on the outside and a faster, magnetically-dominated jet on the 
inside. 

• The era of radiatively-cooled MHD simulations of black-hole accretion disks is beginning. 
This will allow more direct comparison between numerical simulations and observations. 
One desired outcome of this work is an understanding of state transitions in microquasars, 
which may be related to changes in the thermodynamic state of the disk. 
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